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Abstract 



We show that in the loop-erased random walk problem, the exponent 

characterizing probability distribution of areas of erased loops is superuni- 

versal. In (^-dimensions, the probability that the erased loop has an area A 

varies as A~ 2 for large A, independent of d, for 2 < d < 4. We estimate the 

exponents characterizing the distribution of perimeters and areas of erased 

loops in d = 2 and 3 by large-scale Monte Carlo simulations. Our estimate of 

the fractal dimension z in two-dimensions is consistent with the known exact 

value 5/4. In three-dimensions, we get z = 1.6183 ± 0.0004. The exponent 

for the distribution of durations of avalanche in the three-dimensional abelian 

sandpile model is determined from this by using scaling relations. 
64.60.Ak, 05.20.-y, 05.40.+j, 75.10.Hk 
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I. INTRODUCTION 



The loop erased random walk (LERW) problem was first defined by Lawler as a more 
tractable variant of the well-known self-avoiding walk problem. The problem turns out to be 
related to many well-studied problems in statistical physics, but seems to have attracted less 
attention than it deserves. It was shown by Lawler || to be equivalent to a special case of the 
Laplacian self-avoiding walk problem defined by Lyklema et al. ||. Majumdar showed 
that this model is equivalent to the classical graph-theoretical problem of spanning trees on 
graphs, and the g-state Potts model in the limit q — > 0. This equivalence also makes this 
problem related to the abelian sandpile model of self-organized criticality In fact, as we 
shall show in this paper, this model provides a numerically efficient method of determining 
the only unknown critical exponent of the abelian sandpile model in three-dimensions. 

This prompted us to undertake the numerical study of the LERW's in d = 2 and 3 re- 
ported in this paper. We obtain fairly precise estimates of the fractal dimension of LERWs 
in d = 2 and 3. We note the interesting consequence of the scaling theory that the distribu- 
tion of the area of the erased loop has the same exponent 2, independent of the dimension 
d, for 2 < d < 4. In d = 3, the numerical value of fractal dimension of LERWs enables 
us to determine the avalanche durations exponent of the abelian sandpile model, using the 
scaling relations and other exactly known exponents of the model ||. 

A good review of earlier results on the LERW problem can be found in 0. Lawler 
showed that the fractal dimension z of LERWs is 2 for d > 4, and z < (d + 2)/3 for 
(I < 4 0. Recently, it was shown rigorously that in two-dimensions z is strictly larger 
than 1 ||. Using the known exact results about the critical exponents of Potts model from 
conformal field theory, Majumdar was able to prove exactly that z = 5/4 for LERW problem 
in d — 2 [Q, a result which was guessed earlier by Guttmann and Bursill from numerical 
simulations 0. A proof of this result without using conformal field theory has been given 
by Kenyon [IT]]. Using conformal invariance, Duplantier has obtained the exact probabilities 
of no intersection of n LERWs of £ steps starting near each other in two- dimensions, and 
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also the winding angle distribution [[11]]. The distribution of sizes of erased loops was first 



studied in |]T2[ . Priezzhev has used bounds on intersection probability of loop erased walks 
with random walks to show that the upper-critical dimension of the Bak-Tang-Wiesenfeld 
(BTW) sandpile model is 4 0. 

The plan of this paper is as follows: The LERW model is defined in Sec. [J. In Sec. |T|, 
we recall the main points of the scaling theory of the distribution of erased loops ||12|| , and 
apply it to show that the exponent characterizing the distribution of area enclosed by erased 
loop is the same for 2 < d < 4. We determine the behavior of the distribution functions 
for the perimeter and the area of the loop in the scaling limit, for very small or very large 
values of the argument of the scaling functions. The simulation technique and the results 
obtained are described in Sec. [IV]. The exponent characterizing the distribution of durations 
of avalanches in BTW sandpile model in d = 3 is determined in Sec. [Vj, and some concluding 
remarks follow in Sec. [VJ. 



II. DEFINITION OF THE MODEL 



Consider a simple random walk on a <i-dimensional lattice. We start with a particular 
realization W of the random walk having N steps, W = {wq, w\, W2, ■ ■ ., Wjv-ij wn}, where 
Wi is the site reached by the i-th step of walk. We define the LERW £, corresponding to 
W as the path obtained from W by erasing each loop as soon as it is formed. If W has 
no self-intersections, we define C = W. If W has self-intersections, let j be earliest step 
which leads to self-intersection in W, so that j is the least integer such that Wj = Wi for 
some i < j. Then, we obtain a new walk W = {wo, Wi, . . ., Wi, Wj-i, . . ., to/v-i, wn} by 
deleting all steps between i and j, keeping i and deleting j. This process, corresponding to 
loop erasure of the earliest loop formed, is repeated till loops can no longer be found. The 
resulting walk C is the required LERW corresponding to W. This procedure of loop-erasure 
is illustrated in Fig. [1|. 

The length of £ is the number of steps in C We will denote it by n. For a fixed N, n is 

3 



a random variable. We define the critical exponent z of the LERW by the relation that 

(n) ~ N z ' 2 (1) 

for large N, where the angular brackets denote averaging over all random walks of N steps. 
As the root-mean-square end to end distance R is same as for random walks, we have 
R ~ iV 1/2 , and (n) ~ R z . Thus, z is the fractal dimension of the LERW. 

III. SCALING THEORY FOR THE DISTRIBUTION OF LOOP-SIZES 

Let Prob(£, N) denote the probability that a loop of perimeter £ will be erased at the 
JV-th step of the random walk. Let F(£, N) denote the cumulative probability that a loop of 
perimeter £ or greater will be erased at the N-th step of the random walk. We shall study 
the behavior of this function for large N, and write 

Prob(f) = lim Proh(£, N) (2) 

N— »<X> 

and 

F{£) = lim F(£, N) (3) 

N— >oo 

We adopt the convention that if no loop is formed at a step, it will be said to be erasure 
of loop of perimeter 0. With this convention, we clearly have F(0, N) = 1, for all N. 

For d < 4, the mean number of loop-length erased per step tends to 1 for large N. This 
implies that 

oo 

^£Prob(£) = l, ford<4 (4) 
e=o 

For large £, F(£) is expected to vary as a power of £, say as £~ T+1 . However, for a finite 
N, there is a cutoff size £*, and loops of size £ > £* are very unlikely. The cutoff value £* 
varies as a power of N, say £* ~ A^*. This suggests that F(£, N) satisfies the scaling form 

F(£, N) ~ r T+1 f(£/N 4 ') (5) 



The cutoff exponent can be determined by the following simple argument ||12|| : The 
cutoff for the perimeter of erased loops should also vary as (n), the average length of the 
LERW after iV steps. Since this scales as N z ^ 2 [Eq. ([I])], we get <p = z/2. 

The exponent r is also expressible in terms of z. For I -C £*, the total number of loops 
of size > £ for a walk of N steps varies as NF(£), and is much greater than 1. For £ > £*, 
we expect a much stronger decay. For 2 < d < 4, we get a significant number of large loops, 
and thus in this case, we expect that 

NF(jC k ,N)~0(l) (6) 

Putting in the scaling form (|5|), this implies that 

r = 1 + 2/z (7) 

Thus, the scaling form for the distribution of loop perimeters is determined in terms of 
a single exponent z, and is given by 

F{£, N) ~ £- 2 ' z f(£/N z ' 2 ), for £ > 1. (8) 

The scaling function f(x) is assumed to tend to 1 as x tends to zero, and tend to zero 
for large x. We define 

AProb(£, N) = Prob(£, N) - Prob(£) (9) 

If for x near zero, 1 — f(x) varies as x a , we see that keeping £ fixed, and in the limit of 
large N 

AProb(£, N) ~ -K e N~ az/2 (10) 

where Ki is an ^-dependent constant, and the exponent is independent of £. It is easy to 
calculate Prob(2, N) in arbitrary dimension d. The conditional probability of forming a loop 
of perimeter 2 at the iV-th step is 0, if the random walker returned to origin at step (N — 1), 
and it is l/(2d) otherwise (for a d-dimensional hypercubical lattice with coordination number 
2d). Thus 



Prob(2,iV) = ^(l-^) (11) 

where g^ is the probability that the random walker returns to origin after N — 1 steps. In 
(^-dimensions, g n varies as N- d ' 2 for large N. Thus, we see that for large N, AProb(2, N) 
varies as N~ d / 2 . Comparing this with Eq. (|10D, we see that a = d/z. Thus, we get 

f(x) ~ 1 - Kx d/z , for a; near 0. (12) 

[We shall denote an undetermined constant by K. Its value in different equations need not 
be the same.] For other values of I ^ 2, this then implies that 

K e ~ K£ d/z , for t > 1. (13) 

This may be understood as follows: The main deviation of Prob(£, N) from its asymptotic 
value comes from the cases when the LERW at step N—l is of length < £, and the probability 
that walker after N steps is within a sphere of radius i 1 ^ centered at the origin varies as 

£ d/z N -d/2 for £ < N _ 

We can also determine the leading iV-dependence of Prob(£ = 0, N). Since for any non- 
zero £, Prob(£, N) is less than than its limiting value for large N, Prob(0, N) must be larger 
than Prob(0). In fact 

oo 

AProb(0, N) = - AProb(£, N) 

~KN- d / 2 Y,£ [d ~ 2 ~ z)lz (14) 

i 

This summation over i has an upper cutoff proportional to t . In two dimensions, we get 

AProb(£ = 0, N) ~ if (log N)/N (15) 

In three dimensions, the summation diverges as (i*) 1 ^ ~ N 1 ^ 2 . Thus, we get 

AProb(£ = 0, N) ~ K/N, for d = 3. (16) 

For large x, f(x) is expected to decrease as exp(— Kx b ). The exponent b can be de- 
termined as follows: We note that for any constant e < 1, the probability that a loop of 
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perimeter eN is formed at N-th step should vary as exp{— K(e)N} for fixed e and N tending 
to infinity ||15|| . This implies that b — 2/(2 — z), and hence 

f(x) ~ exp(^-Kx 2 ^ 2 ' z) for large x. (17) 

An interesting quantity is the area enclosed by a loop. In two-dimensions, this is straight- 
forward to determine. In three-dimensions, it may be defined as the minimum number of 
plaquettes required to form a simply-connected surface bounded by the loop. In this study, 
we used an alternate, computationally simpler, measure of this area. We simply project the 
loop on to the three coordinate planes, and measure the areas of the projections. If the 
three areas are a 1; 02 and 03, we define the area of the loop to be (a 2 + a 2 , + cl 2 ) 1 ^ 2 - The 
generalization to higher dimensions is obvious. 

Let F(A, N) be the probability that a loop of area greater than or equal to A generated 
at the iV-th step of the random walk. A loop of perimeter £ has a linear size R ~ t x l z and 
an area A ~ R 2 , then, it is easy to see from Eq. (|8]) that for N, A ^> 1 

F(A,N) ~ A~ 1 g(A/N) (18) 

Here also the scaling function g(x) goes to a constant for x — > 0, and decreases rapidly to 
zero for x ^> 1. 

Thus we find the rather unexpected result that the distribution for the area of the loop 
is independent even of z, and hence is the same for all dimensions d, with 2 < d < 4. This 
argument does not work in d = 1, as there Prob(£) decreases exponentially with £, and the 
scaling theory assuming power-law decays fails ||14|| . For d > 4, the LERW behaves as a 



random walk, and for random walks, the area of loop varies as the perimeter of the loop. 
Hence we would expect that the probability that a loop of area A is formed varies as A~ d l 2 
for d > 4. The probability that a loop of area greater than or equal to A is formed, varies 
as A- d l 2+l for d > 4. 

Using the fact that A varies as l 2 l z for t < £*, from Eq. flT2] ) we see that the function 
g(x) determining the finite- N cutoff effects varies as 



g(x) ~ g(0) exp(-Kx d/2 ), for small x. (19) 

For large x, g(x) should vary as exp(— Kx c ), where c is some exponent. Let e be a small 
number <C 1. Using the fact that loops of area eL 2 should decrease only as exp{— K(e)L} 
for fixed e in the limit of large L, we see that c = 1, and 

g{x) ~ exp(— Kx), for large x. (20) 

It is interesting to compare this behavior with that of f(x) for large x. These behaviors are 
consistent only if for a large loop of area A>iV, the average perimeter varies as 

i ~ t{AlNf~ z l 2 (21) 

This behavior should be contrasted with the behavior for A <C N, where the average 
perimeter varies as A z l 2 with no explicit dependence on f. It is interesting to note that 
this scaling law for large loops remains valid even outside the scaling limit for A of order 
N a with 1 < a < 2. For a = 2, it gives perimeter proportional to N, as it should. 

IV. NUMERICAL SIMULATIONS 

The simplest algorithm to simulate the LERW problem on computer is to actually gen- 
erate the trail of a random walk step-by-step on a d- dimensional lattice. At each new step, 
if a loop is formed it is erased. This is straight-forward to implement, but requires large 
memory in large dimensions d, as for simulating a walk of iV steps, one needs to have a 
lattice of linear size N 1 ^ 2 , which means that the required memory increases as 0(N d ^ 2 ). 

The algorithm that makes the most efficient use of memory would store the walk as 
a linked list, keeping only the unerased steps. But there is a memory/CPU tradeoff, and 
the computation time increases as searching for self-intersections is very inefficient in this 
scheme. 

In our simulations we used a hybrid scheme for storing the coordinates of the points 
visited by the walk. We store the coordinates of the LERW in not one, but M lists, where 
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M is a large number. There is a unique hashing rule which assigns a site to one of the 
lists, so that checking for intersection has to be done only within one list. To see if a point 
already belongs to the LERW, we have to search only in the list corresponding to the point 
in question. The best choice of M is 0(£*), as then each list has 0(1) entries. With this, 
we were able to simulate a two-dimensional walk of 2 33 steps in about 3 hours 16 minutes 
using ~ 60 Mb of memory on a 350 MHz Pentium-II machine. In three-dimensions, a walk 
of 2 29 steps took about 24 minutes and ~ 300 Mb of memory on a similar machine. 

Simulations were carried out for total walk length N of 2 r steps, with r = 25, . . . , 29 
for two-dimensional walks and r = 24, . . . , 28 for three-dimensional walks. To eliminate the 
initial transients, we collected the statistics of loops only after discarding the first N/2 steps. 
In addition, ensemble average was taken over 10 3 distinct realizations of random walks in 
each case. For the two-dimensional case we also simulated a small number of walks for r up 
to 34. 



A. Two-dimensional LERW 



In Fig. H| we show the plot for the cumulative distribution function F(£,N) JT^]. We 
plot £ 2 ' Z F(£, N) versus £. There is a significant deviation from simple power-law behavior 
for very small £, and for large £. For I < £*, the data fits well to the functional form 
given by Eq. (|T2|). In the small £ regime, the leading correction is a correction to scaling. 
Incorporating this, we fit the data to the form 

Fo(i,N) = -^exp [-C^/N)-"' 2 ] (l + ^ (22) 

where C*2 is related to the fractal dimension via C*2 = 2/z. 

The best fit values of all the parameters in Eq. ( p2"D are tabulated in Table |. We note 
that C5 is 1 within our error bars. Furthermore, the exact value of C*2 is also known to be 
8/5. As a result, one more set of values were estimated for the parameters by constraining 
C*2 and C5 to these values. The parameter values thus obtained are also tabulated in Table |. 



The fit is rather good for all t > 10. Statistical fluctuations are large for £ > 10 4 , as there 
are not many such loops generated. 

In Fig. |H we have plotted the AF(A, N) versus A for different values of N, and also shown 
the best fit using the fitting form Eq. (|22|) with i substituted by A. While estimating the 
parameters we constrained C2 and C5 to 1. This is because the exact value of C2 is known 
to be 1 and unconstrained value of C5 turns out to be 1 within error bars. This allows a 
better estimate of the remaining parameters. The estimated best fit values of parameters 
for this data set are tabulated in Table |. It is clearly seen from Fig. |3| that the scaling form 
fits the data very well in nearly the entire range. 

We obtained more accurate estimates of Prob(£, N) for i < 100 by taking an ensemble 
average over 10 9 different realizations of the random walk. In Fig. |], we have plotted the 
variation of iVAProb(£, N) versus log(iV) in two dimensions for i = 0, 2, 4 and 6. We see 
clearly that while iVAProb(0, N) has a linear variation with log(iV), for other values of £, 
this tends to a limiting constant value for £ ^> 1. 

B. Three-dimensional LERW 

The distribution of loop-sizes for the three-dimensional walks by perimeter is shown in 
Fig. p. The format of presentation is exactly the same as in the previous subsection. We fit 
the data to the form given by Eq. (^). From the figure it is seen that this scaling form fits 
the entire data very well for i > 10. The best fit values of parameters in this equation are 
tabulated in Table |. We find that in this case the best-fit value of the correction to scaling 
exponent C5 turns out to be 0.86, clearly different from 1. 

Our estimate of the best fit value of the fractal dimension z gives 

z = 1.6183 ±0.0004, d = 3 (23) 

This value is not very sensitive to the choice of parameters C%, C3, C4, and C5. The error bar 
on z gives our subjective estimate of errors of extrapolation. This should be compared with 
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the value z = 1.623 ± 0.011 obtained by Guttmann and Bursill |J. Because of the larger 
fractal dimension of walks, for the same value of N, there are significantly more longer loops 
generated in d = 3 than in d — 2. As a result, we see power-law scaling over roughly 5 
decades of £ in Fig. [5] compared to that of about 4 decades of £ in Fig. 0. 

In Fig. H we show the plot for F(A, N) the cumulative distribution function for loop 
area. We have plotted AF(A, N) versus A for different values of N. An unbiased estimate 
of C2 from the best-fit gives a value 1.00000 ± 0.00005. Thus, we put C2 to be exactly 1 
and estimated the remaining parameters by fitting the scaling form given by Eq. fl22|) with 
£ substituted by A. From the figure it is clearly seen that this form approximates the entire 
data very well for A > 10. The estimated values of parameters is tabulated in Table |. Here 
also, the exponent in the correction to scaling term turn out to be different from 1. 

In Fig. |7| we have replotted the data of Fig. § with AF(A, N) plotted against (A/N) 3 / 2 . 
We see that the curves are approximately linear for small A/N, verifying the theoretical 
prediction of Eq. [19]. For larger values of A/N, the slope decreases as expected from Eq. |2(J 

In Fig. we have plotted the variation of AProb(£, N) versus N for £ = 0,2, and 4. 
The data was obtained by averaging over 10 9 different realizations of 100-stepped walks. For 
Prob(£), we used the values from the N = 2 28 simulation. We see good agreement with the 
predicted 1/N variation for £ = 0, and 1/N 3 ^ 2 variation for £ = 2 and 4. 



V. RELATION TO EXPONENTS OF THE SANDPILE MODEL 



The sandpile model of Bak-Tang-Wiesenfeld is defined as follows [fL7|] : We consider a 
hypercubical lattice of linear size L in d dimensions. At each site is a non-negative integer 
which gives the "height" of the pile at that point. The system is driven by adding a grain 
of sand at a randomly chosen site, thereby increasing the height of pile at that site by one. 
If the height at any site exceeds (2d — 1), it topples, and its height decreases by 2d, and one 
grain is transferred to each of its neighbors. If this makes some other sites unstable, they 
are toppled in turn, until all sites are stable, and then a new grain of sand is added. 
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Bak et al. observed that in the steady state of such a pile, adding a grain gives rise to a 
sequence of topplings, and the size of such avalanches is a random variable with a power-law 
tail. Determining the exact values of the exponents characterizing these tails has been the 
main theoretical problem in the area of self-organized criticality. 

The model in d = 1 is rather trivial, and does not show simple power-law tails of 
avalanche sizes, as most avalanches are large fl8|]. In d = 2, unlike in d = 1, most avalanches 



are finite, but they involve multiple topplings of sites. A theoretical understanding of this 
case remains incomplete ||19|| . For d > 4, mean-field description of avalanche propagation 



is adequate, and the corresponding exponents are the same as of sizes of clusters in critical 
percolation theory [f[3] . 

The BTW model for d = 3 does not suffer from the problems caused by multiple top- 
plings. It is thus the simplest undirected model for studying self-organized criticality with 
nontrivial (non-mean-field) critical behavior. In this case, multiple topplings at a site occur 
with very low probability, and the avalanche clusters are found to be compact, with fractal 
dimension 3. Then, simple scaling arguments |2t]|J show that if the probability that there 
are exactly s topplings in an avalanche in a system of linear size L is Prob(s | L), which 
satisfies the simple finite-size scaling form 

Prob(s | L) ~ S - a h{s/L b ) (24) 

then we must have a = 4/3 and 6 = 3. 

The theoretical assumptions that go into the scaling argument have been checked exten- 
sively in simulations, but a rigorous theoretical proof is not yet available. Since the number 
of distinct toppled sites is assumed to be proportional to the number of topplings, we see 
that the probability that an avalanche has Sd distinct sites toppled also varies as s^ 4//3 . 

As a check on the scaling theory, note that the probability that avalanche reaches a 
distance R scales as the probability that number of topplings is greater than R 3 , hence 
as 1/R, which also agrees with the known result about expected number of topplings at a 
distance R. 

12 



The only exponent which this simple argument does not give is the exponent for the 
duration of avalanches. But the propagation of avalanches occurs along spanning trees path 
in the equivalence between the sandpile model and spanning trees 0. Hence, the duration 
T of an avalanche must vary with its linear extent as T ~ R z . And the z is for spanning 
trees, which is the same as the z we used for LERWs. The knowledge of z, thus, allows us 
to estimate the exponent for duration of avalanches: the probability that the duration of 
avalanche is greater than T varies as T~ v , where 

y = i/ z = 0.61795 ±0.00015 (25) 

VI. CONCLUDING REMARKS 

We have already noted that the LERW problem is very suited for numerical studies. In 
the two-dimensional case, we have collected data for over 10 3 realizations of walks with N 
up to 2 29 . Thus the numerically determined loop-size distribution is an average over more 
than 8 x 10 10 loops (only about 31.25% of the steps taken form non-trivial loops on square 
lattice). For the three-dimensional case, the corresponding number is 2.8 x 10 10 loops (only 
about 21.17% steps form non-trivial loops on cubic lattice). The quantity which corresponds 
closest to loop-erasures are avalanches in the sandpile model (more correctly, subavalanches) 
PHI . Clearly, simulation of the Abelian sandpile model with equal number of avalanches is 
not possible with available computing machines. 

Secondly, our simulations are done on an effectively infinite lattice, and there are no 
complicated boundary effects to complicate the analysis of data. Corrections due to finite 
size of system show up only in the finiteness of the number of steps N of the random 
walk. This seems to be well described by simple finite-size scaling theory. If we wanted 
to determine the exponent z using the sandpile model, or the spanning trees, the largest 
system sizes accessible would be much smaller. 

The dimension-independence of the exponent characterizing the distribution of areas of 
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erased loops for 2 < d < 4 is rather unexpected. The exponent does depend on dimension 
for d > 4. We have been unable to find a more transparent proof of this result. 

We would like to thank S. N. Majumdar for his critical reading of an earlier version of 
this paper. 
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TABLES 

TABLE I. Estimated values of various parameters corresponding to Eq. ( |22| ) for variation of 
loop perimeter and area for two- and three-dimensional LERWs. The values without error bars 
are fixed during estimation of other parameters. 







Ci 


c 2 


c 3 


c 4 


c 5 


2-D 


£ 


0.3533 ± 0.0004 


1.5997 ± 0.0005 


1.1 ±0.1 


1.58 ±0.03 


1.000 ± 0.007 




£ 


0.35385 ± 0.00025 


8/5 


1.1 ±0.1 


1.56 ±0.03 


1 




A 


0.127316 ±0.000015 


1 


9.8 ±0.7 


0.494 ± 0.003 


1 


3-D 


£ 


0.1527 ±0.0003 


1.2359 ± 0.0003 


2.8 ±0.3 


1.69 ±0.02 


0.86 ± 0.05 




A 


0.1312 ±0.0002 


1 


35 ± 7 


0.142 ±0.005 


0.394 ± 0.015 
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FIG. 1. An illustrative example of the loop-erasing procedure: The random walk 
a-b-c-d-e-- • — p starts at a, and ends at p. The erased loops are shown by thin lines. Note that 
at the point o, while the random walk path intersects itself, the LERW encounters no intersection, 
as the loop f-g-h-i-f has already been erased. 
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FIG. 2. Variation of £ 2 ^ Z F(£, N) with I for two-dimensional LERW. Solid lines represent the 
best fit of Eq. ( p2|) with parameter values (including estimate of z) given in Table |. 
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FIG. 4. Variation of iVAProb(£, iV) with N for £ = 0, 2, 4 and 6 in two dimensions. For better 
visibility, the data points for I = 2, 4, and 6 have been displaced vertically by 0.69, 0.59, and 0.57 
respectively. 
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FIG. 5. Variation of £ 2 ^ Z F(£, N) with £ for three-dimensional LERW. Solid lines represent the 
best fit of Eq. ( p2|) with parameter values given in Table |. 
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FIG. 7. Variation of AF(A,N) with 100(A/N) 3/2 for threedimensional LERW showing the 
behavior of the exponential correction term. Thick straight line represents the best fit of Eq. (|2 
with parameter values given in Table | 
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